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Abstract. In the Bondi-Sachs formulation of General Relativity space-time is 
foliated via a family of null cones. If these null cones arc defined such that their 
vertices are traced by a regular world-line then the metric tensor has to obey 
regularity conditions at the vertices. We explore these regularity conditions when 
the world line is a time-like geodesic. In particular, we solve the Einstein equations 
for the Bondi-Sachs metric near the vertices for axially symmetric vacuum space- 
times. The metric is calculated up to third order corrections with respect to a 
flat metric along the time-like geodesic, as this is the lowest order where non- 
linear coupling of the metric coefficients occurs. We also determine the boundary 
conditions of the metric to arbitrary order of these corrections when a linearized 
and axially symmetric vacuum space-time is assumed. In both cases we find that 
(i) the initial data on the null cone must have a very rigid angular structure for the 
vertex to be a regular point, and (ii) the initial data are determined by functions 
depending only on the time of a geodesic observer tracing the vertex. The latter 
functions can be prescribed freely, but if the vertex is assumed to be regular they 
must be finite and have finite derivatives along the geodesic. 



PACS numbers: 04.20.-q, 04.20.Ex, 04.20Jb, 04.25D-, 04.25. Nx 
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1. Introduction 

The pioneering work of Bondi and Sachs led to a rigorous understanding of 
gravitational waves at large distances from a compact source (Bondi et al. 1962, 
Sachs 1962). They used a coordinate chart (a; , x 1 , x 2 , x 3 ) = (it, r, 6, 4>) for the 
metric that was adapted to out-going null cones and where only one coordinate, the 
areal distance r, varies along the out-going null geodesies generating these cones. 
The other coordinates are the 'retarded time' u, and the standard spherical angles 
9 and <f>. Assuming a Minkowskian observer at infinitely far distances from an 
isolated gravitating object they showed that such a system can lose mass (energy) 
only via gravitational radiation. Their results were confirmed by Newman and Unti 
(1962) who used the Newman-Penrose formalism (Newman and Penrose 1962, 1963) 
to demonstrate the advantage of null foliations in General Relativity in order to 
characterize gravitational waves. 

An observer at null infinity can only keep track of what is happening in the future 
or past of a source that emits gravitational waves. However, if one is interested in 
studying the emission process itself, e.g. by means of numerical simulations, one needs 
to consider conditions at the source. Tamburino and Winicour (1966) were the first to 
impose the boundary conditions for the gravitational field at a finite distance on the 
null cone, i.e. they placed a Minkowskian observer on a time-like world-tube of finite 
space-like radius, e.g. the radius of a star. This world-tube null formalism, where out- 
going null geodesies are attached to a world-tube, forms the basis of many numerical 
codes that are currently used to solve the Einstein equations in a domain extending 
from a finite-size world-tube to null infinity (see Winicour 2012a, for a recent review). 

When the radius of the world-tube becomes zero, the world-tube degenerates 
into a world-line, i.e. into the time-like curve of an inertial observer tracing both the 
vertices of out-going null cones and the origin of the coordinate system. This so-called 
Fermi observer is a non-rotating observer in a rectangular Minkowskian coordinate 
system along the world-line. Isaacson et al. (1981) integrated the Einstein equations 
for the Bondi metric (Bondi 1960) when the world-tube had zero radius. The boundary 
conditions for the Bondi metric were chosen such that it approaches flat space values 
at the origin. To assure this behavior they restricted the fall-off of the metric functions 
towards r — to certain positive powers of r. The work of Isaacson et al. motivated 
Gomez et al. (1994) and Siebel et al. (2002) to solve numerically the vacuum Einstein 
and the Einstein-fluid equations for the Bondi metric in axial symmetry from the 
vertex to null infinity. In their integration schemes, the boundary conditions for the 
Bondi metric included the lowest order curvature contribution. Both used a correct, 
but ad hoc, ansatz for the boundary conditions near the vertex, which they did not 
provide for all variables, however. 

In this work we systematically investigate the boundary conditions of the Bondi 
Sachs metric near the vertex. The vertex of a null cone is the focal point where 
all null geodesies generating the cone either converge to or emanate from. From 
the mathematical point of view, the null cone is not differentiable at its vertex, 
and consequently derivatives of the metric and therefore the curvature tensor cannot 
be calculate there. In particular, any tensor tensor that is expressed in terms of 
coordinates, like the Bondi-Sachs coordinates, adapted to the null cone cannot be 
expanded in terms of a Taylor series with respect to these coordinates at the vertex. 
Therefore, the boundary conditions for any fields expressed in Bondi-Sachs coordinates 
arc a priori not known at the vertex. In principle these boundary conditions can 



The Bondi-Sachs metric at the vertex of a null cone 



3 



be chosen freely. However, the vertex is the origin of an inertial observer along 
the geodesic, who implies regular boundary conditions. The assumption of such an 
observer at the vertex also requires that the vertex is a regular point for the considered 
tensor fields on the null cone (Dautcourt, 1967). A regular point & (the vertex) for a 
function / (e.g. a metric component) on a open domain W (i.e. the null cone without 
the vertex) is defined as a point on the boundary , where / has a Taylor expansion 
based on & into °U (Freitag and Busam, 2005). If the vertex is a regular point in the 
null cone, we will call it a regular vertex. 

To find the proper boundary conditions of fields on a null cone at its vertex, we 
have to make additional assumptions: 

(i) Assuming that the vertices of the null cones are traced by a world-line of a Fermi 
observer, we require this world line to be a time-like geodesic since then the 
acceleration of the curve vanishes, which greatly simplifies the analysis. This 
geodesic should be a regular time-like geodesic, i.e. its tangent vector should 
nowhere vanish along the curve. 

(ii) The geodesic should be a contained in a convex normal neighborhood (Hawking 
and Ellis, 1973). Such a neighborhood, T say, has the property that two distinct 
points in T can be connected by a unique geodesic and at any point, in T 
Riemann normal coordinates, z a , can be established such that the metric in the 
neighborhood of & can be written as 

.gab=»?ab-^ acbd (^)z C Z d + ... , (1) 

where r\ a p — diag(— 1, +1, +1, +1) and 7£ ac bd(^) is the Riemann tensor 
evaluated at In this regard, a metric is called a regular metric at a point 
^ G T, if it is invertible at & and its determinant is finite and nonzero at 

These two assumption are crucial to find the boundary conditions of the Bondi-Sachs 
metric at the vertex of the null cone. These two assumptions are necessary to find the 
boundary conditions. Assuming a regular metric of an inertial observer moving along 
a time-like geodesic, we will use coordinate transformations to find a Bondi-Sachs 
metric at the time-like geodesic. This Bondi-Sachs metric will not be regular at the 
vertex because the angular base vectors are not defined there. Nevertheless, by this 
procedure we will obtain the desired boundary conditions for the Bondi-Sachs metric 
corresponding to an observer in an inertial frame at the vertex, because the initial 
metric was a regular metric. To this end we will address four questions: 

(I) What are the generic requirements of the boundary conditions of the Bondi-Sachs 
metric near the vertex of the null cone if the vertex is a regular point in the null 
cone? 

(II) How can these boundary conditions be determined? 

(III) What are the explicit boundary conditions for the Bondi-Sachs metric of an 
axially symmetric vacuum space-time? 

(IV) What do the regular boundary conditions imply for the initial data on a null cone 
with a regular vertex for axially symmetric vacuum space-times? 

To find the answers to these questions, we construct a Fermi normal coordinate 
system (Misner and Manasse, 1963) along the time-like geodesic. Based on this 
rectangular coordinate system, we then calculate a metric on a null cone where the 
radial coordinate is an affine parameter. This null metric is subsequently transformed 
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to a Bondi-Sachs metric, by changing the radial coordinate over to an areal distance 
coordinate. This has the advantage that the Einstein equation split into a hierarchy of 
hypersurface equations and evolution equations (Bondi et al. , 1962, Winicour, 2012a) 
which simplifies their numerical treatment on a null cone. If an affine parameter is 
used instead, one of the hypersurface equations contains an additional time-derivative 
which is numerically more challenging (Winicour, 2012b). Alternatively, we could 
start with a regular metric along the time-like geodesic in deDonder coordinates as in 
Thorne (1980), Thorne and Hartle (1985), Zhang (1986) and calculating the boundary 
conditions for a Bondi-Sachs metric via an appropriate coordinate transformation. 
Since these authors derived their metric expansions for vacuum space-times, only, we 
do not uses their expansions to answer question (I), since assuming a vacuum space- 
time restricts the generic properties of the boundary conditions (see Sec. 3.1). Second 
(and more importantly) starting from a Fermi normal system is simpler because the 
Fermi normal coordinates are geometrically adapted to the properties of the convex 
normal neighborhood. We will later compare our results with those of Zhang (1986), 
who obtains an expansion of the metric in the 3 + 1 formulation of space-time. 

The behavior of arbitrary fields near the vertices of null cones was first discussed 
by Penrose (1963), who realized that these fields should be expanded in terms of 
normal coordinates defined at the point of origin of the expansion. His argument 
was formal providing only a sketch of the expansion. Friedrich (1986) considered the 
expansion for the case of the conformal vacuum Einstein equations, and, like Penrose, 
used the spin frame formalism to derive the connection and the curvature variables. 
As we are interested in the conditions at the vertex of the null cone in the Bondi-Sachs 
formulation, where the curvature is given as a function of the metric, the method of 
Penrose (1963) and Friedrich (1986) does not directly apply. 

Ellis et al. (1985) discussed regularity conditions at the vertices of past null cones 
which are traced by a time-like geodesic of a cosmological (comoving) observer. Their 
work motivated Poisson and collaborators (Poisson, 2004, 2005, Poisson et al. 2006a, 
2006b, 2010, 2011) to use light cone coordinates at a world-line to obtain the metric of 
a non-rotating black hole moving along the world line. As both Ellis et al. and Poisson 
and Vlasov used an affine parameter instead of an areal distance along along the null 
rays, they would have had to transform the former into the latter by an appropriate 
coordinate transformation to obtain the Bondi-Sachs metric. 

In principle, we could have used the results of Ellis et al. and Poisson et al. as 
a starting point to answer questions (I)- (IV), but our aim is to present a complete 
and pedagogical derivation of the boundary conditions of a Bondi-Sachs metric at 
the vertex from a regular metric of an inertial observer along a time-like geodesic. 
With this in mind, we follow the approach of Ellis et al. (1985, App. A), but we 
provide additional pedagogical steps in the derivation and generalization of their 
results. Regarding Poisson and collaborators, we do not need the machinery of a 
bi-tensor theory or the Synge world function (Synge, 1960) to find a null metric with 
an affine parameter as radial coordinate. However, we will recover the results of Ellis 
et al. and Poisson et al. in our calculations. 

Choquct-Bruhat et al. (2010, 2011) and Chrusiel and Jezerski (2010) investigated 
the Cauchy problem of General Relativity for initial data given on the null cone, 
and Chrusciel and Paetz (2012) took these studies up comparing the results with 
those of Friedrich (1986). They analyzed the 'vertex problem' raising questions about 
uniqueness and existence of solutions. Although these results are mathematically 
very interesting, they are too general to help in formulating boundary conditions 
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for a Bondi-Sachs metric, and in answering question (III) and (IV) above. We note 
that regularity issues at the origin of a curvi-linear coordinate system are also of 
importance in the 3+1 formulation of General Relativity. The first in-depth analysis 
of this problem for axially symmetric systems is due to Bardeen and Piran (1983). 

The article is organized as follows: After defining our notation in the reminder of 
the introduction, we derive in Section 2 the limiting behavior of a Bondi-Sachs metric 
from a regular Fermi normal coordinate system, and answer question (I) concerning 
the generic properties of the Bondi-Sachs metric. In Section 3, we discuss and respond 
to question (II)-(III). In Section 4, we conclude with a summary and discussion of the 
results, and we respond to question (IV). 



1.1. Notation 

We assume the existence of a smooth, four-dimensional manifold M A with a Lorentzian 
metric g and signature (— , +, +, +). The flat metric is always denoted with 77. Greek 
indices (a, j3, ...) and the small Latin indices (a, b, ...,f) run from to 3, while the 
small Latin indices (i,j,k,...) can assume the values 1, 2, or 3, respectively. Capital 
Latin indices (A, B, C ...) label angular directions on a sphere and take the values 
2 or 3 for the coordinates (6, 4>) or (y = —cos9, <fi), respectively. The Kronecker 
symbol is Sij — diag(+, +, +). Unless stated otherwise, the Einstein summation 
convention holds. We use five distinct symbols for the four-dimensional space-time 
coordinates: (i) y a are arbitrary coordinates; (ii) z a are Riemann normal coordinates; 
(iii) y a are Fermi normal coordinates; (iv) x a are affine null coordinates; and (v) x a arc 
Bondi-Sachs coordinates. Partial derivatives are denoted by a comma, and covariant 
derivatives by the V— symbol. For the Christoffel symbols r^„, the Riemann lZ a 'p 
and Ricci tensor TZ^, we follow the convention of Misner et al. (1972). Moreover, 
we use geometrized units where the speed of light and the gravitational constant are 
equal to unity. Associated Legendre polynomials of the first kind, P[ n (y), are defined 
as in Jackson (1999), i.e. 



2. A Bondi Sachs metric derived from a regular metric 

2.1. A regular metric in Fermi normal coordinates 

Let (^ 4 , g) be the four-dimensional space-time. Suppose (^ 4 , g) contains a simply 
convex normal neighborhood V with a regular time-like geodesic c(r), where c" is 
its unit, time-like tangent vector and r the proper time. Since Y is a convex normal 
neighborhood and c(r) is contained in T, a regular metric g a p can be found along 
c(t). Let y a be a set of four-dimensional arbitrary coordinates along c(r) such that 
the metric g a p is regular and can be expanded along c(r) like 

=g a f) , Xv") +9ap,Ki , y 1 + jjffa/wn _y K1 y K2 + -, (2) 

c(t) Z 



c(r) 



c(t 



where the coefficients are evaluated along c(t). Because of arguments given in section 
2.3 below, it is sufficient to consider the expansion only up to the quadratic term. 

At any point on c(r) we choose an orthonormal tetrad e a (r) which is parallel 
propagated along c(r), i.e. V c ( r )e a = 0. The time-like base vector of the coordinate 
system e is tangent to c(r), i.e. e (r) = d/dr. On any point of the geodesic c(r), we 
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send out space-like geodesies b(r, n, £) which are parametrized by an affine parameter 
£ and point into the direction n = n'e^r), i.e. n° = and n = d/d£\ c i T y The affine 
parameter £ of the geodesies b(r, n,£) is defined to be zero along c(r). This specifies 
the coordinates y a as Fermi normal coordinates y a (Misner and Manasse, 1963, Ni 
and Li, 1979), which are given by 

y° := r , y [ := £n [ , £{y { ) := , (3) 

where aV/dr = and n 1 = y 1 /£ are the direction cosines of the space-like geodesies 
b{T,n,t). Note that the three-dimensional coordinates y l are Riemann normal 
coordinates % for every value of t on c(r). Constructing the coordinates in that 
way assures that the metric along c(r) is the Minkowski metric, because the base 
vectors are orthonormal along c(r). In addition, it implies that the first order partial 
derivatives, g Q/ 3, Kl , vanish along c(t), because of the parallel transport equation of 
n 1 along c(r) and the geodesic equations of c(r) and 6(r, n, £). The particular form 
of the second partial derivatives at the geodesic c(r) was first derived by Misner and 
Manasse (1963), and the corresponding metric reads 

3oo(y a ) = - 1 - ^oioj , W)H + <W] 3 ) , (4a) 

c(t) 

3ok(y a ) = - l^oikj f (r)yW+0([y a } 3 ) , (46) 

6 c(t) 

5km (2/ a ) - <5ij - Jftikjm , (r)2/V + G([y a ] 3 ) , (4c) 

6 c(t) 

where 7^ a bcd are the Riemann normal components of a Fermi observer along the time- 
like geodesic c(r). These components are space-time invariants, i.e. knowing all 
twenty independent Riemann normal components along c(t) allows one to construct 
uniquely (up to spatial rotations) a metric up to quadratic terms of a power series 
expansion with respect to Fermi normal coordinates along the time-like geodesic in a 
sufficiently small neighborhood T of the geodesic. 



2.2. An affine null metric 

The null coordinate x° := f w labels null cones whose vertices are along the time- 
like geodesic c(r), where f w is equal to the proper time r along the geodesic, i.e. 
Tw\c(t) — T - For points not on c(r), f w is constant along parametrized null geodesies 
that emanate from c(t). We introduce the three spatial coordinates (x 1 , x 2 , x 3 ). The 
radial coordinate x 1 := £ is a positive affine parameter for null rays emanating from 
c(t), where £ = 0. The two coordinates x A := (x 2 , x 3 ), which are constant both along 
the time-like geodesic c(r) and the null geodesies, label the direction angles of the 
null rays. In the following, we will use the notation when we refer to an arbitrary 
point on c(r) that is also the vertex of an arbitrary null cone. For £ to be an affine 
parameter, it has to obey the condition 

k a VJ. = l , (5) 

where V a is the covariant derivative with respect to the metric in Fermi normal 
coordinates. Equation (5) holds, if we define £ as in (3) and choose the null vector 
k a (y h ) ■= (1, y l /£) in the Fermi normal coordinate system. Then £ does not only 

| For a thorough discussion of Riemann normal coordinates, see e.g. Schouten (1954), Thomas 
(1991), and Iliev (2006). 
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measure the space-like distance from the geodesic c(r) in a Fermi frame, but it also 
gives the positive affine distance between a vertex 6" and points on the respective null 
cone of 6. We note that an affine parameter along a curve is generally defined only up 
to the transformation A — a£ + b, where a and 6 are arbitrary constants. However, our 
choice £ — on c(r) implies b = 0. From a = sign(a)|a| and the definition of £ (3) one 
sees that multiplying the Fermi normal coordinates y 1 by \a\ corresponds to scaling 
these coordinates, i.e. without loss of generality we can set \a\ = 1. Hence, null rays 
fc(A) emanating from c(r) are most generally parametrized with the affine parameter 
A = w£, where w = ±1. We choose w = 1 for future-pointing null rays along c(r) and 
w = — 1 for past-pointing ones by imposing the normalization condition 

lim c a fc a = -w (6) 

along c(t), where c a is the tangent vector of c(r). 

A null geodesic fc(A) in the null-cone f w = const emanating in x A — direction, 
admits an expansion in Fermi normal coordinates y a with respect to the affine 
parameter A of the form (using A = w£ and w 2 = 1 ) 



£ 



2 



y a (^) = f w 8\ + w£k a y b (x A ) + -F a y b (f w , x A ) + — G a y b (f w , x A ) +0(£ 4 ), 



wt 



3 



3! 



(7) 

where the coefficient functions are evaluated along c(r). Setting £ = in (7) shows 
that it describes the time- like geodesic c(r) with the tangent vector c a = S a T , 
and that the tangent vector of the null geodesies k a = (1 / 'w)dy a / 'd£\ c i T ^ does not 
depend on f w , because the x A are constant along c(r). If c(r) were no geodesic, 
an additional dependence of f w through the directional vector k a must be taken into 
account (Newman and Posadas, 1969). Hereafter, we use the null vector k & in the form 
fc a = (l, ^'(x' 4 )), where n 1 is a three-dimensional unit vector being parametrized by 
the angles x A . 

The geodesic equations of the null rays k(w,£) read in Fermi normal coordinates 
d 2 V a _ ( d] dy h dy c 

-^--T hc (y ) — — , (8) 
and further differentiation with respect to £ gives 
<*V F a , g ,dy b dy c dy d 

Inserting (7) into (8) and (9), and equating the thus obtained relations along c(r) by 
setting £ = leads to 

F a (? w ,x A )=0 , 

G a (f w ,x A )= -Tg cd (f w )k b (x A )k c (x A )k d (x A ) . 

c(t) 

In Fermi normal coordinates G a vanishes, because the derivative rg c d can be expressed 
as a linear function of the Riemann tensor 7^ bcd , which is antisymmetric in the last two 
indices, and these two indices are contracted symmetrically. Hence, the parametric 
representation of the null rays k(w,£) that emanate from the time-like geodesic c(r) 
reads in the Fermi frame 

y a (x a ) = f w 5% + w£k h {x A ) + 0(£ 4 ) . 
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The metric in affine null coordinates, x a , is given by the coordinate transformation 

9a ^)= gah [y^)]^{xn%{xn , 

and the Jacobian (dy a /dx a ) between y a and x a by 
flu 11 P 2 w£ 3 

^{^)=5\ + -F% w (f w ,x A ) + —G a fw {f w ,x A ) + 0{£ i ) , (10a) 

^V) = wk a (x A )+(D(£ 4 ) , (106) 
r)n a f 2 in? 3 

J^(^) = w£k a A (x B ) + yF* A (f w , x B ) + —G a A {f w , x B )+0(^) , (10c) 

which contains the derivatives of the coefficients F a and G a with respect to the null 
coordinates. Although these coefficients vanish at £ = 0, their derivatives are not 
necessarily zero there, i.e. we have to calculate the derivatives first, and subsequently 
evaluate them in the limit I = 0. Using the first derivatives of the Christoffcl symbols 
in a Fermi normal coordinate system as in Ni and Li (1979), we calculated the metric 
9ap{x a ) in affine null coordinates near the time- like geodesic c(r): 

g 00 (f w , £, x A ) = - 1 - £ 2 A(f w , x A ) + 0(£ 3 ) , (11a) 
g 0A (f w , £, x A ) = - ^w£ 3 B A (f w , x A ) + 0(£ 4 ) , (116) 
9i„(f w , £, x A ) = -w5\ , (11c) 

g AB (fw, I, x A ) = £ 2 q AB (x c ) - jS AB (f w , x A ) + 0(£ 5 ) , (lid) 

where q A B{x c ) = 8ijn 1 A r\} B is a unit sphere metric around G, because Tn}(x A ) is a unit 
vector labeling points on a sphere, and the coefficients 



A(f w ,x c ) =^ 0a0 b (f w )k a (x u )k b (x°) , (12a) 



— T^-OaOb 


c(r) 




= "^-Oabc 


c(r) 




— 7^-abcd 


cfr) 



(? w )k a {x c )k h A (x c )k c (x c ) , (126) 



S AB (f w ,x°)=K ahcd (f w )k a (x c )k b (A (x c )k c (x c )k d B) (x°) (12c) 

c(t) 

depend on the Riemann normal components 7?. a bcd , Sjw) evaluated along the 
geodesic c(r), i.e. 7^ a bcd| , Af w )=TZ a b c d(T w ). We point out that the components <?i M 
are exact, and as they are equal to a constant, they do not require an order symbol. 
In particular, |gio| = 1, because £ is an affine parameter. 

The components of the affine null metric (Ha)- (Hd) obtained by us agree with 
those in the literature (Ellis et at 1985, Poisson et al. 2010). While Ellis et al. (1985) 
calculated a null metric at a past null cone only, we derived the metric for points 
both in a past and future null cone, which is achieved by setting w = 1 (future) or 
w = — 1 (past) in our solution. The difference between our work and that of Poisson 
et al. (2010) concerns the computational approach of how to find the affine null metric 
along c(t). Poisson et al. (2010) use the Synge world function (Synge, 1960) to obtain 
the metric, whereas we use an explicit coordinate transformation based on a regular 
Fermi normal coordinate system. Both approaches are legitimate and it is a matter 
of taste which one prefers. 
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2.3. The Bondi-Sachs metric 

Using the affinc null metric g^x 01 ) derived in the previous subsection, we calculate 
a Bondi-Sachs metric and the corresponding Bondi Sachs coordinates are labeled by 
x a . The null coordinate x° = t w and the angular coordinates x are defined like 
their barred counterparts, f w and x A , respectively. The radial coordinate x 1 := r is a 
positive areal distance coordinate, i.e. surfaces dr w = const and dr = const have the 
area Airr 2 . This also requires that the determinant of the 2- metric g A B divided by r 4 
does not vary with time and radius. 

The areal distance r is defined by (Jordan et al. (1961), Sachs, 1961, Newman 
and Penrose, 1963) 

fc M V M r = r6(fc) , (13) 

where Q(k) := iV^fc^ 1 is the expansion rate of null rays with tangent vector k^. 
Inserting the null vector k^(x a ) — wS^^ into equation (13) we obtain 

r(xn =l(l- jS(f w , x A ) + 0(£ 3 )) 1/4 , S(f w , x A ) := q AE ' (x A )S AB (f w , x A ) . (14) 

The equation shows that the affine parameter is equal to the areal distance r both for 
\l\ s=y and in a flat space-time. The latter holds, because for a flat metric S vanishes 
due to its dependence on the Riemann tensor. 

Wc define a coordinate transformation between the affine null coordinates x^ and 
the Bondi-Sachs coordinates x^ according to 

r w =f w , r(x^=l-^S(? w ,x A ) + 0(e 4 ) , x A = x A (15) 

Applying this coordinate transformation to the metric components (Ha) - (11 d) results 
in the following Bondi-Sachs metric near the geodesic c(r): 

gaaixP) = - 1 - r 2 A(r w , x c ) + 0(r 3 ) , (16a) 
9oa(x*) = - ^wr 3 B A (r w , x c ) + G(r 4 ) , (166) 

2 

gi^n = -w[l+ r -S( Tw , x c ) + O(r 3 )]s tl , (16c) 

gABix") = r 2 q AB (x c ) - r — [s ab {t Wi x c ) - ^q AB (x c )S(T w , x c )] + 0(r 5 ) . (16rf) 

Note that the above metric is indeed of the Bondi-Sachs type, because the determinant 
of gAB /t 2 is a function of x A , only. 

Traditionally, the Bondi-Sachs metric is written in the form (Bondi 1960, Sachs, 
1962) § 



ds 2 



-c^+^drl - 2we 2f3 dT w dr + r 2 h AB (dx A - U A drJ) (dx B - U B dr w ^ , (17) 

where the determinant of the 2-metric, dct(h A B), is a function of x A only. Because 
of this restriction, the tensor h A s possesses only two degrees of freedom. Using the 
standard spherical coordinates as angular coordinates, x A — (9, </>), allows us to write 
h A B in the form (van der Burg, 1966) 

h AB dx A dx B = e 27 cosh(2(5) d9 2 + 2 sin 6 smh(25)d9d(f) + cosh(2(5) sin 2 6 d<j) 2 . (18) 

§ (1) Bondi and Sachs, and also most numerical relativists (see Winicour (2012) for a review) use 
retarded time u := r+i as time coordinate; (2) for computational convenience, we renamed the 
original metric function V/r as cxp(2$ + 2/3) in the metric component goo{x a ). 
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By linearizing Hab with respect to 7 and S, it can be seen that 7 and 6 correspond 
to the two components of a two-dimensional transverse and trace-less tensor, e.g. the 
tensor xab m Chusciel et al. (1998, 2002). According to (17) the Bondi-Sachs metric 
5a / s(a ;M ) is not defined at r = 0, because the four-dimensional volume element \J~ r g 
vanishes there. 

For convenience we introduce the quantity & <G {7, S, (3, $, U A } to abbreviate 
the set of the six functions describing the Bondi-Sachs metric. To find the behavior 
of & at the vertex, we assume that & can formally be expanded into a power series 
at r = 0, i.e. 

„2 



JT(^) = &M(t w , x a ) + r&V>(T w , x A ) + ^&W(t w , x a ) + 



2! 



(19) 



where the coefficient functions ^^{t w , x A ) are evaluated along the geodesic c(r). 
Inserting the expansion of each element of & into (17) and (18) results in a series 
expansion for every component of the Bondi-Sachs metric in terms of the elements of 
& . A comparison of these series expansions with the expression given in equations 
(16a) - (16d) yields the following relations between & and the contracted Riemann 
normal components A, Ba and Sab'- 



7 (x a ) 
5(x a ) 
f3(x a ) 



See{r w , x A ) sin 2 6 - S^{t w , x a ) 



12 

r 2 S, 



sin 



0(r 3 ) 



G 



+ G(r 3 ) 

Sill u 

S(r w ,x A ) + 0(r 3 ) , 



(20 a) 
(206) 
(20c) 



U A {x a ) = 



AC 



$(x Q ) = — A(t, 



B c (r w ,x D ) + 0(r 2 ) 



1 



-S(t w ,x d )\ +0(r 3 ) 



(20 d) 
(20 e) 



Equations (20a)-(20e) show how the traditionally-used Bondi-Sachs metric functions 
& behave to lowest order near the vertex of the null cone, when the vertex coincides 
with the origin of a Fermi normal coordinate system along c(r). Hereafter, these 
lowest order corrections are also referred to as first-order corrections, ^i, of the 
Bondi-Sachs metric with respect to the flat metric at the vertex along the time- 
like geodesic or briefly first-order corrections of the Bondi-Sachs metric. These first 
order corrections ^\ arise at 0(r) for U A and at 0(r 2 ) for /?, 7, S, and where as 
higher-order correction coefficients ^„ are expected at 0(r n ) for U A and at 0(r n+1 ) 
for (3, 7, <5, and $, respectively. 

Since equations (20a)-(20e) were derived form a regular metric at the vertex and 
the vertices were assumed to be on a timeline geodesic, we deduce (without solving the 
Einstein equations) the following regularity requirements on a power series expansion 
of the traditionally-used functions which assure that the Bondi-Sachs metric is 
regular at its origin: 

(i) The power series of the metric functions & in r must start at r = with a certain 
positive power of r, i.e. 7, 5, (3, $ are of 0(r 2 ) and U A of 0(r). 

(ii) The angles x must parameterize topological spheres centered at the vertex and 
the radial coefficients of & must show a specific angular behavior determined by 
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the contractions of the Riemann normal components || with the tangent c a , the 
null vector fc Q , and k a A along c(r). 

In Sec. 3.2.1, we will show that the topological requirement is crucial for 
regularity, since it is also possible to have non-regular solutions at the vertex, 
when this assumption is not imposed. 

(Hi) The radial expansion of & have specific numerical factors in the radial expansion 
coefficients, for example the factor 1/12 in (20a). 

Two additional important properties of the functions & at the vertex cannot 
be inferred from the first-order corrections in (16a)-(16d), because they manifest 
themselves only at second- or third-order deviations from the flat metric at r — 0. 
While the first-order corrections originate from the second derivative of the Fermi 
metric at the origin through a coordinate transformation, the second and third-order 
corrections result from the third and fourth derivative of the Fermi metric. To learn 
more about the two missing properties we consider a qualitative argument, which 
also avoid tedious calculations of the corresponding coordinate transformations. We 
recall that in Fermi normal coordinates the Christoffcl symbols vanish along c(r) , and 
that in any coordinate system, where the Christoffel symbols vanish, the following 
correspondence can be made between third and fourth partial derivatives of the metric 
and the Riemann tensor (Schouten, 1954): 



where d n g is the n^-partial derivative of the metric, L(-) a linear functional of 
its arguments, (Riem) the Riemann normal components, and V™ the n t/l -covariant 
derivative, respectively. The explicit dependence between the third and fourth-order 
partial derivatives of the metric and the Riemann tensor in Fermi normal coordinates 
is given in Ni and Li (1979), Dolgov et al. (1983), and Ishii et al. (2005). 

To calculate the Bondi-Sachs metric including the second and third-order 
corrections, the corresponding afline null metric must be determined first. Following 
and extending schematically the procedure of section 2.2 to higher order in the 
approximation reveals that the afline null metric involves time derivatives of the 
Riemann normal components at the vertices on c(r). These derivatives arise because 
covariant derivatives of the Riemann tensor in the Jacobian are contracted, like e.g. 
k^VaA = fc a fc b fc c V a 7^obOc with the null vector k a , which is given by k a (x A ) = 
5 a Tw + n l (x A )S a { . Starting at second-order corrections, these time derivatives show 
up in hierarchical manner, i.e. the second-order corrections contain first-order time 
derivatives of first-order corrections to the flat metric at the vertex, like e.g. A, Tw . 
Similarly, the third-order corrections contain second-order time derivatives of first- 
order corrections, like e.g. A, TwTw . Moreover, there are first-order time derivatives of 
second-order corrections in the third-order correction coefficients, like e.g. the time 
derivative of the second-order correction n 1 (x A )\?iA. 

At this stage, it is instructive to compare where time derivatives of the 
Riemann normal components occur in the metric expansion using different coordinate 
conditions. In table 1, we list the occurrence of the contracted Riemann normal 

| Physically speaking, this requirement means that if the space-time (^f 4 , g) is curved, the path of 
a null ray emanating from the vertex is affected in the neighborhood of the vertex (i.e. the Fermi 
observer) by the curvature (the Riemann tensor and its covariant derivatives) at the vertex. 



d 3 



g = L V(Riem) 



(21a) 



d 4 g = L V 2 (Riem), (Riem) 2 



(216) 
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Table 1. Occurrence of the contracted Ricmann normal component ■S'rjiOj" 1 "'' , 
their contracted spatial covariant derivatives and their respective time derivative 
along the time-like geodesic in the metric component goo m affinc null coordinates, 
Fermi normal coordinates and de Donder coordinates. The order of the correc- 
tion coefficient V n corresponds to the power of r n+1 in the radial expansion of goo • 



order of 


affine 


Fermi 


deDondcr 


correction 


null 


normal 




coefficient 


coordinates 


coordinates 


coordinates 


% 












r R- i0 j ; k n'n i n k 


r R-n )j0 . ik n t nin k 










^3 




7liQjQ. k in % riin k n l 


TZiOjO; k in l nin k n l 




d{TZiQj . k n % nin k )/dT w 















components IZoiOj & n d their spatial covariant derivatives in the first three correction 
coefficients of goo to a flat metric on the geodesic for an expansion in Fermi-normal 
coordinates, in deDonder coordinates and in affine null coordinates, respectively. The 
behavior for the Fermi normal coordinates can be taken from Ni and Li (1979), 
Dolgov et al. (1983), or Ishii et al. (2005), and the one for the deDonder coordinates 
from Zhang (1986, eq. 3.26). In Fermi normal coordinates there are no time 
derivatives of the Riemann normal components or their spatial derivative, whereas 
in deDonder coordinates the lowest (second) order time derivative occurs in the 
third order correction coefficient. In affine null coordinates, however, we recognize 
a hierarchical order of the time derivative of the Riemann normal components and 
their spatial covariant derivatives^. 

Following the same line of arguments as for the time derivatives, one can show that 
the third-order correction coefficients depend on the square of the Riemann normal 
components. This behavior of the affine null metric near the vertex was also observed 
by Poisson and Vlasov (2010), who derived an affine null metric for vacuum space- 
times. Since the Bondi-Sachs metric and the affine null metric are related by a 
transformation of the radial coordinate only, everything that we said above about 
the affine null metric applies to the functions J^, too. 

As normal coordinates exist along c(r) and as a power series of a metric in 
normal coordinates can be transformed into a power series of a metric in Bondi-Sachs 
coordinates, regularity implies the following restrictions on the power series expansion 
coefficients of & at the vertex: 

(iv) The correction coefficient (n > 1) depends on the time derivatives of order 
(n — k) of the correction coefficient Cfe, (1 < k < n) and these time derivatives 
must be finite. 

% In deDonder coordinate, these time derivatives also occur in hierarchical order, i.e. the ^ 2 fc 
correction coefficients of the metric along the geodesic depends on the (2k) th time derivative of 
the k th covariant derivative of IZoiOj ■ 
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(v) Non-linear coupling occurs at lowest order near the vertex in the third-order 
correction coefficient of & , i.e. at 0(r 4 ) for 7, 5, (3, $, and at C(r 3 ) for U A . 

The requirements (i) — (v) state the general properties of the boundary conditions for 
the Bondi-Sachs metric functions & , when the vertex is a regular point in the null 
cone which is traced by a time-like geodesic. As these properties have been derived 
without the solution of the Einstein Equations (they can be considered as "kinematical 
conditions" ) , they depend on the solution of Einstein equations (which is a "dynamical 
condition"). 

We summarize the general requirements on the boundary conditions of the 
functions & at & in the following Vertex Lemma: 

Let the Bondi-Sachs metric functions, & ', be represented by power series expansions 
with respect to the areal distance coordinate r at the vertex of a null cone on a time-like 
geodesic. If the coefficients of this radial expansion obey 

(1) the general properties (i) — (v), 

(2) and the Einstein equations, 

then this power series expansion of & can be used to formulate boundary conditions 
for the metric functions & at the vertex of the null cone, and the vertex is a regular 
point in the null cone. 

In the next section, we explicitely calculate the boundary conditions for & at 
the vertex for axisymmetric vacuum space-times. We choose these space-times, 
to complete and extend the boundary conditions as used by Gomez et al. (1994) 
and Siebcl et al. (2002), and to recover and justify their ad hoc assumptions 
on the boundary condition employed in their numerical algorithms. In addition, 
axisymmetric space-times are the simplest + ones to show all the relevant addressed in 
the vertex lemma. 

3. Solutions for axially symmetric space-times 

A four-dimensional axially symmetric space-time (^# 4 , g) is a space-time containing 
a time-like 2-surface srf consisting of points that are invariant under the action of a 
one-parametric, cyclic group G that is isomorphic to SO(2) (Carter, 1970, Stephani 
et al. 2003, and references therein). In particular, since the metric is invariant under 
the action of G, i.e. there exists a Killing vector field £(</>) in (^# 4 , g), where <j> is 
the parameter of the group G. The vector field £(<fi) is constant along the orbits of 
the group action, and the Lie derivative JP^^g^ vanishes along the curves generated 
by £(</>)■ It can be shown (Carot, 2000, and references therein) that the 2-surface 
si ', the so-called axis of symmetry, is auto-parallel. Since (^# 4 , g) is endowed with 
a metric, the fact that srf is auto-parallel is equivalent to srf being totally geodesic 
(Spivak, 1999). As srf is time- like, the axis of symmetry contains a family of time- 
like geodesies. Since these geodesies distinguish themselves from other geodesies in 

due to their axial symmetry, we call them axial geodesies. They are the natural 
choice to trace the origin (vertex) of a Bondi-Sachs coordinate system in an axially 
symmetric space-time. 

Hereafter, let c(r) be an axial geodesic that is normalized in such a way that 
r is the proper time. Given an orthonormal tetrad e M (r) along c(r), we define a 

+ Spherically symmetric space-times are of no use here, because they possess no angular structure^, e. 
property (ii) of the vertex lemma cannot be demonstrated. 
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Fermi normal coordinate system y a along c(r), where d/dr is tangent to c(r) and r is 
the time coordinate. We further choose three mutually orthogonal space-like vectors 
d/dy l at every point on c(r), which are parallel transported along the axial geodesic. 
Of these three vectors, the two vectors d/dy 1 and d/dy 2 are normal to the time-like 
2-surface srf , while the vector d/dy 3 is tangent to it. When y° = const and y 3 = const, 
the coordinates y 1 and y 2 label points of Killing orbits of G in .J? A whose fixed-points 
are located on the axial geodesic c(r). When instead y° — const and y 1 = = y 2 
holds, the coordinate y 3 labels points on the symmetry axis srf ' . In the Fermi normal 
coordinate system introduced above the Killing vector has the form (Carot, 2000) 

e{y°) = y 2 S\ - yH\ . (22) 

In the following, we will present two approaches to determine the behavior of 
the Bondi-Sachs metric functions & near c(r). In the first approach, in Section 3.1, 
we solve the Einstein equations in the Fermi normal coordinate system to obtain the 
Ricmann normal components, which arc then used to calculate the Bondi-Sachs metric 
functions from the expressions (12a) - (12c). In the second approach, in Section. 3.2, 
we expand the metric functions & into a power series obeying both axial symmetry 
and the limiting behavior of & near r = as given in (20a) - (20e). Subsequently, we 
solve the vacuum Einstein equations in Bondi-Sachs coordinates for the power series 
coefficients. 



3.1. Lowest order non-tivial boundary conditions for the Bondi-Sachs metric derived 
directly from the Fermi metric 

The vacuum Einstein equations are 72^ = 0. In axial symmetry, every tensor T has 
to obey the Killing condition, i.e. the Lie derivative Jz?£(0)T has to vanish. If we 
apply this condition to the Riemann normal components 72 a b c d in a Fermi normal 
coordinate system using the Killing vector (22), we find the following relations among 
the non-zero components of IZabcd'- 

72-0101 = 72-0202 ,72-0113 = 72o223 ,72o3l2 = — 272oi23 ,721313 = 722323 , 72o303 ,72i212- 

(23) 

For the further discussion we conveniently combine these non-zero components into a 
set I e {A, B, C, D, E, F} with 



A := 72oioi , B := 72o303 , C :— 72oii3 , D := 72oi23 
The Ricci tensor reads in Fermi normal coordinates 
/ 



E := 72i2i2 , F := 72i 3 i 3 . (24) 



72 a/3 (y a ) 



-2A-B 



-2C 







A-UE + F) 






\{E- 




F) 



-2C 



B-F 



+0(y a ) ,(25) 



c(r) 



where the functions I are evaluated along the axial geodesic c(r). We note that the 
function D\ , . is not determined by the Ricci tensor. Hence, it is not determined 

I c(t) j ' 

by the vacuum Einstein equations in Fermi normal coordinates, and thus can be 
prescribed freely at the vertices of the null cones. Solving the vacuum Einstein 
equations for the zeroth-order coefficient of the expansion (25) results in the solution 



B 



c(r) 



= -2A 



C(T) 



c\ 



c(r) 



= , E\ 



C(T) 



= 4A 



C(T) 



Ic(t) 



= -2A 



c(r) 



(26) 
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which depends only on the function yll , freely specifiable along c(r). 

Next we set A\ c ^(t w ) := 672 (t w ) and D\ c ^(t w ) := 6S 2 (t w ), and calculate the 
Bondi-Sachs metric up to the first-order correction to the fiat metric for an axially 
symmetric vacuum space-time using (12a) — (12c) and (20a) — (20e). We find 

l{x a ) = r 2 j 2 (T w )P 2 (y) + 0(r 3 ) , (27a) 

f3(x a ) = (D(r 3 ) , (276) 

U e (x a ) = - 4n2(T w )P}(y) + 0(r 2 ) , (27c) 

$(x a ) = -6n 2 (T w )P°(y) + 0(r 3 ) , (27d) 

and 

S(x a ) = r 2 5 2 (T w )P 2 (y) + 0(r 3 ) , (27e) 

U+{x a ) = - 4rS 2 (r w ) (^4) + 0(r 2 ) , (27./) 



sin( 

where y := — cos 9 and -P; m (j/) are the associated Legendre polynomials of the first 
kind. The functions 72 and 5 2 determine, respectively, the electric and magnetic part 
of the Weyl tensor at the vertex. Moreover, they also define the lowest order curvature 
contributions in the norm and twist of the Killing vector w.r.t. their flat space values 
at the vertex, respectively. 

Equations (27a)-(27/) give the lowest non-trivial order boundary conditions for 
& at a regular vertex for axially symmetric vacuum space-times. Equation (276) 
shows that (3 does not behave as 0(r 2 ) as expected by property (i) of the vertex 
lemma. This also demonstrates that if we had started with the expansions of Zhang 
(1986) to find the generic properties of the Bondi-Sachs metric at the vertex, the most 
general limiting behavior of j3 would be wrong. In fact it can be shown that for axially 
symmetric perfect fluid space-times (3 is of 0(r 2 ) at the vertex. 

In principle, it is possible to calculate also the next order correction coefficients 
with the above approach, but this involves contractions of the covariant derivatives 
of the Riemann normal components over five and more indices with the null 
vector fc a and its angular derivatives, e.g. the 0(r 3 ) coefficient of 7 depends on 
(V a 7?-bcdo) k a k h k c A k d k c B . Thus, we describe in section 3.2 another approach that 
avoids the tedious calculation of these contractions, and also is easily extendable to 
higher order. 

3.2. Regular boundary conditions for the Bondi-Sachs metric derived from a power 
series of ^ 

To find axially symmetric solutions for the Bondi-Sachs metric (17), we transform the 
Killing vector (22) to Bondi-Sachs coordinates x a , i.e. £ a (x a ) = 5°^. The Killing 
equations J^^u\g a ^(x^) — imply that the Bondi-Sachs metric functions & do not 
depend on the coordinate <j>. To facilitate the computations, we employ the coordinate 
transformation y = — cos6> in the metric (17). Points on the axis of symmetry, i.e. the 
poles, are given by y = ±1. We further define the auxiliary function s(y) := \J\ — y 2 , 
and anticipating from (20 a) - (20 e) we assume the following expansions for the metric 
functions 

N 

l(T w ,r,y) - ^( W r)"+ 1 7„ +1 (r tu ,y)+o[( U ;r) Ar + 2 ] , (28a) 

n=l 
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6(T w ,r,y) =^{wr) n+1 5 n+1 {T w ,y) + 

n=l 

JV 

f3(r w ,r,y) = Y.( wr ) n+1 ^+i{T w ,y) + 

n=l 

JV 

n=l 
JV 

U A (r w ,r,y) = J2( wr ) nU n(r w ,y) + 

n=l 

where the set of expansion coefficients ^ n := {7^+1, S n +i, (3 n +i, &n+i, U A } is 
calculated at r = 0, N is the order up to which the field. Note, this expansion only 
respects property (i) of the vertex lemma, and we do not impose other restrictions at 
this stage. After calculating the vacuum Einstein equations for the coefficients <£ n , we 
end up with coupled partial differential equations for the ^„ with respect to u and 
y. We solve these equations in general, and then restrict their solution to be regular, 
i.e. finite, at the boundaries. This will provide us with the boundary conditions for 
the Bondi-Sachs metric at a regular vertex . These regular boundary conditions at 
6 will then also comply with the general properties stated in the vertex lemma. In 
the following, we refer to ^ n as the n th — order correction of the Bondi-Sachs metric 
functions & with respect to a flat metric at 6 or simply the n th — order correction of 

In Section 3.2.1 and 3.2.2, we solve the vacuum Einstein equations, which can 
be grouped in Bondi-Sachs coordinates into six so-called main equations, three 
supplementary equations, and one trivial equation, respectively (Bondi et al. 1962, 
Sachs, 1962, van der Burg, 1966, Winicour, 2012). The six main equation can be 
split further into two evolution equations for the transverse-traceless part of the 2- 
metric /iab(7, 5), and four hyper-surface equations for the variables (3, U A , and $. 
Furthermore, from the twice contracted Bianchi identities follows the lemma (Bondi 
et al. 1962, Sachs, 1962, Tamburino and Winicour, 1966): If the main equations hold 
on one null cone and if the optical expansion rate of the null rays does not vanish 
on this cone then the trivial equation is fulfilled algebraically and the supplementary 
equations hold if they are fulfilled at one radius r. Consequently, we need to consider 
only the main equation to find the solution of the Einstein equations at the vertex &. 
The supplementary equations then provide a check of this solution. The Ricci tensor 
components and the quantities derived from it that appear in the supplementary, 
hyper-surface and evolution equations determining the corrections ^ n are given in 
Appendix A . 



(wr) 



JV+2 



(286) 



(wr) 



JV+2 



(28c) 



(wr) 



JV+2 



(28 d) 



(wr) 



JV+l 



(28e) 



3.2.1. Boundary conditions depending on 'ifi, ^2 o,nd ^3 According to property (v) 
of the vertex lemma and the expansion of a metric in normal coordinates, see e.g. 
relation (216), we expect the non-linear coupling of the ^-coefficients to happen at 
lowest order in the ^3-coefficients. Indeed, using (A. 7), lZ rr = and omitting terms 
0(r 3 ) gives 

/3 2 = , (29a) 
& = , (296) 
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1 r 



(72) 2 + (5 2 f 



(29c) 



Equation (29c) shows that the ^3 correction of (3 depends quadratically on the solution 
of the ^1 correction. 

To find the solutions for the n G {1,2,3}, coefficients we solve first for ^\ 
and ^2 setting the constant and linear coefficients of the r— series in (A. 7) - (A. 12) 
equal to zero. We next insert the solution of the ^-correction into the quadratic 
coefficients of the r— series of the Ricci tensor components in (A. 7) - (A. 12) and solve 
for the ^-coefficient. 

Inserting the solution (3 2 = = f3 3 into (A. 8) - (A. 11) and utilizing TZ a p =0wc 

find 



= 2Uf + 2(s 2 72 ) ty + [5C/| + 3( S 2 73 ),y] rw' 1 + 0(r 2 ) , 

= 2Uf + 2s- 2 (s 2 5 2 ), y + \hUt + Ss-V^Jrw- 1 + 0(r 2 ) , 
= - 12$ 2 + W? <y + 2 S - 2 {s i l2 , y ) y - 4 72 



+ 



wr + 0{r 2 ) 



(30 a) 
(306) 

(30 c) 



- 16$ 3 + 6Ul y + 2 S - 2 (s 4 l3 ,y). y - 4 73 j 
= - 127 2 + 3s 2 (s- 2 C/f), 2/ + [l2 72iTro -24 l3 + 4s 2 (s- 2 Uy),y]wr + 0(r 2 ) , (30(f) 
= 12<5 2 - 'is 2 ul y + ( - 125 2 ^ + 245 3 - As 2 ul^wr + 0(r 2 ) . (30c) 
Combining algebraically the zeroth-order term in the r— expansion of (30a) and (30(f) 



gives 



72 



= s 2 72, w - 2J/72,,, + (6 
while the same procedure 

= s 2 73,j/j/ - 2J/73,!/ + (l2 - -^73 - 5 72 , 



(31) 
(32) 



while the same procedure yields for the first-order terms 

4 

7 2 

and leads to similar equations as (31) and (32) when applied to (306) and (30e), 72 
being replaced by S 2 , and 73 by £3, respectively. Hence solving (31) and (32) does 
not only give 72 and 73, but also provides the structure of a solution of S 2 and 63. 
Inserting the ansatz 7 2(t, ui , y) = ■y 2 (T w )S(y) in (31) yields 

4 



= 72(t to ) 



2 \ d 2 S n dS 
dy dy 



6 



S 



(33) 



l-j/2, 

Since 72 is an arbitrary function, the term in the parenthesis has to vanish, which also 
provides us with an associated Legendre differential equation for S (y) of the form 



d 2 S 



dS 



S , 



(34) 



dy 2 dy L 1 — y 2 - 

with I — 2 and m = 2. This equation holds for the Legendre polynomials of the first 
kind, P" l {y) 1 and of the second kind, Q™(y), respectively. Hence the most general 
solution of (34) involves two integration constants Si and S 2 , i.e. 

S(y) = S 1 Pr(y) + S 2 Qr(y) . 

Since the associated Legendre polynomials of the second kind, Q^{y), are singular at 
the poles y = ±1, we set S 2 = 0. Hereafter, we will always set integration constants 
connected to Legendre polynomials of the second kind equal to zero. Thereby we 
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guarantee regular solutions at the poles and point (ii) of the vertex lemma, since the 
P/™(y) are well defined for all angles ye [— 1, 1]. A regular solution of (31) assuming 
separation of variables is 

72 Cr™, y) =l2(T w )P%{y) , (35) 

where we absorbed the integration constant S\ into the arbitrary function 72 (t w ). To 
solve (32), we make the ansatz 



k=2 



(36) 



where the coefficients 73 (t w ) depend on the time at the vertex and are determined by 
the differential equation (32). Inserting (36) into (32), reveals that ^(t^) := ji(r w ) 
is freely specifiable function, j^(t w ) — 5/6(drf2/dT w ), and all remaining coefficients 
73 (t w ) must be zero. Therefore, the solution for 73 is 



7s(w) = ~f 3 (T w )Pi(y) + I \^{t w )\ Pi(y) • 

lUT w J 



(37) 



As already mentioned above, we can derive similar partial differential equations as 
(31) and (32) for <5 2 and S3, i.e. we may assume a solution for 82 and S3 as 

S2(r w ,y) = S 2 {t w )P 2 2 {y) , 



S 3 (T w ,y) = S 3 (T w )P$(y) + 



dS2 
dr,, 



(r w ) Pi{y) , 



(38) 
(39) 



\}s{y) 



where both 52(j w ) and S 3 (t w ) are arbitrary functions of t w . 

To solve for U± and we insert the solutions for 72, 73, <5 2 , and S3 into (30 a) 
and (306) and solve in the constant and linear term of the r— series for the respective 
variable. This yields 

U?(t w , V ) = -4y(r w )P2 1 (y)s(y) , 

Uf(r w ,y)= -4S2(r w )P^y)s- 1 (y) , 

The solutions for 3> 2 and $3 are found inserting the solution for 72, 73, Uf, and 
U% into Eq. (30c) and solving the constant and linear term of the r— series for $ 2 an d 
$3, respectively. This gives 

$2(T w ,y) = -6j 2 (T w )P^(y) , 
$3(r w ,y) = -l2%( Tw )P°(y)- 



U y 2 {T w ,y) = -{&h(T w )P£(y)+2 
U$(r w ,y)= -{6S3(r w )P 3 1 (y) + 2 



(40) 
(41) 

(42) 
(43) 



UT7/1 



(44) 
(45) 



Having obtained the complete solutions for the coefficients ^1 in (29a), (35), (38), 
(40), (41), and (44), and %f 2 in (296), (37), (39), (42), (43) and (45), we now determine 
the ^-coefficients of the metric function We insert the solutions for 72 and 82 
into Eq. (29 c), which gives 

If ~ ~i r „ -1 2 

Pa{t w , y) = ^\[l2(T w )] 2 + [S 2 (t u 
= {[12(t w )] 2 + [5 2 (t w )] 



}[pi(y) 



\pliy) - Y b Pl(y) 



(46) 
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This equation shows that the solution can be expressed again as a linear combination 
of associated Legendre polynomials. We insert the remaining & 2 coefficients into (A. 8) 
- (A. 12), and then set these equations equal to zero. Considering only the quadratic 
coefficients results in the following equations: 



= 9U$ + 4( S 2 74 )^ + 324jys 4 [( 72 ) 2 + (6 2 f 
= 9U$ + 4s- 2 (s 2 S i ), y , 

= - 20$ 4 + 2s- 2 (s 4 7ity ) - 4 74 + Wl v - 9(1(% 4 - 32y 2 + 12)( 72 ) 2 



+ 72(10?/- lly 2 + l)((5 2 



240y S 2 73 , Tl „ + 540s 4 



2 / — 2 

= 40(5 4 - 5s 2 Ul y - AQs z 8 2 ,r m r w - 2AQys z 8 z . Tw + 720s 2 (5 272 . 
Combining (47a) with (47 d), and (476) with (47e) yields 

4 



=- 40 74 + 5s 2 (s- 2 Ul) y + 40 72 , T „ Tli 
- 720s 2 7 2 , 

„2 ; 



(72) 2 + {52? 



(47a) 
(476) 

(47c) 

{Aid) 
{Ale) 



= s 2 -u,yy - 2 yii . y + (20 - ^) 74 - 18s 2 { 72 , TiuTro + 6j/ 73 , r „ - 9[( 72 ) 2 - (<5 2 ) 2 ]}(48) 

= s 2 S it yy - 2yS i!y + (20 - ^jSi - 18s 2 (s 2 . T ^ + 6y V» - 18 72 <5 2 ) . (49) 

Since the square of an associated Legendre polynomial can be expressed by a linear 
combination of Legendre polynomials with the same degree of m but varying I, we 
assume that both 7 4 and 04 obey an expansion in terms of Pf{y) as 

00 00 
lA{T w ,y)^^l{{T w )P^{y) , 8 A {T w ,y) ^^5{{T w )Pf{y) . 

1=2 1=2 

Inserting this ansatz into (48) and (49) provides us with the solutions for the 
coefficients 74 and <5 4 , respectively, 

74(W) = l4{r w )Pi{y) + ^ [j^{t w )\ P!{y) 



+ I 



d 2 l2 

dr 2 



{l~w ) 



h{T w ,y) = 5 4 {T w )P 2 {y) + 
' 3[d 2 S 2 



9 

10 



Y {[l2{r w )} 2 -[S2{r w )} 2 }jP 2 { y ) 

{r w )]Pi{y) 



(50) 



^3 
dr,,, 



7 



dr?, 



{t u 



54~ ~ ] 
- —12{t w )5 2 {t w ) >P 2 {y) , 



(51) 



where "fi{r w ) :— jf{r w ) and 5 4 {t w ) := Sf{T w ) are freely specifiable functions. While 
(47a), (476), and (47c) provide the solutions for and <E> 4 , respectively: 

U y 3 {r w ,y) = ( [8 74 (r w ) + ^ 7 2 (r) + g^r)] P}{y) + 4 \^{r w )] P£(y) 



.dr u 



+ {2ict^)-y[^-)] 2 }- 2 ^}^ 



(52) 



Ut{r w ,y) 



-{ 



SS^Pl^+Al^^Jp^y) 
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+ 



16 d 2 6 2 , , 96^ ~ i w 
21 ~d^( Tw > ~ —^(t w )S 2 (t w )^ P 2 (y) 



20 7 4(t„) + ^[ 72 (t„ 



216 r ~ , 



4 c? 2 72 
7^" 



24 r 
(r w )- y {[72( 



)] 2 

2 + fe)] 2 } 



s-\y) , 



^73 

dr..- 



f Fft(r.)] a -f[fc(r.)] ! 



(53) 

^ 3 °(y) 

(54) 



Equations (46), (50) - (54) constitute a solution for the coefficients of 
The solution for the corrections ^2 and ^3 demonstrates points (m), (iv) and (u) 
of the vertex lemma, as the expansion coefficients contain specific numerical factors 
(point (in)). Time derivatives (point (iv)) of the lower order coefficients ^\ and ^2 
occur hierachically ordered, i. e. ^2 depends on the first time derivative of 'rfi whereas 
^3 depends on the second time derivative of ffi and on the first time derivative of 
^2- The non- linear coupling (point (v)) in the correction coefficients occurs at lowest 
order in ^3, since ^3 depends on (^f) 2 . If one uses the solution for ff n ,n — (1,2,3) as 
boundary conditions for the Bondi-Sachs metric in numerical simulations for axially 
symmetric space, the free functions 72 and S 2 must be at least twice differcntiable with 
finite derivatives, 73 and 63 must be differcntiable with finite derivatives, whereas 74 
and 64 must be continuous. Thus, requiring the vertex to be a regular point in the 
null cone rigidly fixes the boundary conditions fbrj;he metric at the vertex. The only 
freedom left is the choice of the functions 7„ and S n , n£ {1,2,3}. 



3.2.2. Boundary conditions that are linear in ^„ and valid for n > 1 To investigate 
the boundary conditions of & further, we consider (28a) - (28e) for an arbitrary value 
of N. Since the metric functions & are power series in r with coefficients c £ n , the 
corresponding Ricci tensor can be written as 



N 



y) = J2r n K$(T w ,y) 



i.e. as a power series in r, too. As the coefficients 1zffi will, in general, not depend 
linearly on we linearize them with respect to %? ra to firstly simplify the calculations, 
and secondly to find the qualitative dependence of the hierarchical order of the time 
derivatives addressed in property (v) of the vertex lemma. The linearized Ricci tensor 
components for the main equations of the vacuum Einstein equations are given in 
Appendix A. 2. Introducing the operator 

~l(l ■ 



U(f) ■■= (i-y 2 )f, y 



2yL 



1 - 



f 



(55) 



for a function / that depends on y, allows us to write some of the upcoming equations 
in a more compact way. In particular, if f(y) — P 2 (y), we find 



L, 



P?(y) =b\P](y) 



b\ :=l(l + l)-j(j + l) 



(56) 



where there is no summation performed over j on the right hand side of the first 
expression. Note that the associated Legendre polynomial Pf(y) commutes * with 
the operator Lj, because b l l = 0. 

* A similar relation as (56) could be derived for the Legendre polynomials of the second kind, Qf(y), 
but we do not investigate this further here as these polynomials arc not defined at y = ±1. 
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From (A. 13) and the vacuum Einstein equation izffl = 0, we deduce 
p n+1 {T w , r,y) = , Vn> 1 . (57) 
Inserting the solution for f3 n into (A. 14) - (A. 20) and combining = from (A. 14) 
and (A.15) with TZ^ = and TZ^ = from (A.17) and (A.18), respectively, yields 
= L 2 ( 72 ) , = L 2 (<5 2 ) . (58) 

Since the associated Legendre polynomial P$(y) commutes with L 2 , regular 
solutions of (58) for ye [—1, 1] may be assumed to have the form 

72(t™, y) = ll{T w )Pi(y) , S 2 (t w , y) = S 2 2 (r w )P^y) . (59) 
where ^) 2 (t w ) and S 2 (t w ) are arbitrary functions depending on the time t w along the 
geodesic c(t). Given the solutions for 7 2 and <5 2 , and from TZfy — 0, we derive 
U?(T w ,y) = -4j 2 2 (T w )P}(y) S (y) , U? (t w , y) = -^(t^P* (y^ 1 (y) . (60) 

Setting T^2D) m (A.16) equal to zero and inserting 7 2 and Uf leads to 
^2(T Wl y)^-^fl{T w )P^(y) , (61) 
which shows that the solution for the lowest order correction ^1 obtained with the 
linear approach is equivalent to that of the generic approach in the previous section. 

To find the solutions for n > 1, we assume hereafter n > 1. Combining 
TZ 1 ^ = from Eqs. (A.14) and (A.15) with = and = from 

Eqs. (A. 19) and (A. 20), respectively, gives rise to two equations that determine the 
solutions for 7„+i and <5„+i: 

= L n+ i(7„+i) + a„+i7n,r ro , = L n+ i(5 n+ i) + a n+ i5 ntTw , (62) 
where 

._ 2(n + 2)(«-l) 

(n + 1) ' ( ' 

The expressions in (62) show that the solutions for 7„ and S n obey the same type of 
equation, i.e. knowing one of the solutions gives us the form of the other solution, 
too. 

We treat the solutions for 7„ and S n together by defining l£ := Jk^ A y + 3kO~ A ,p 
with k > 2, whereby (62) becomes 

= L„ +1 (7^ +1 )+a„ +1 7^ . (64) 
Making the ansatz 

n+l 

In + i(r w ,y) = E J n+i.k{r w )P 2 k {y) , (65) 

k=2 

which we insert into (64) and obtain for n = 2, 3, and 4 

= h 3 (I A ) + a 3 I 2 A Tm = b\l A 3 Pi + [b\l A 2 + a 3 d -^j Pi , (66a) 

= U(l£) + a,I A Tw = b\l A A Pt + (b\l A 3 + a 4 ^) Pi + (b 2 J A 2 + «4^) P| , 

(666) 

l5 7A n2 , (tA 7A , „ dI 4A \ D 2 , f,3 7 A , „ dI 4.3 



= L 5 (It) + a 5 l£ Tw = b\If*P£ + I b\l£ A + a 5 -^\Pi + \b\l£ 3 + a 5 ^ I P 



+ \b\iL + a^)pi mo 
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Since the associated Legendre polynomials depend on the arbitrary angle y and 
arc non-zero in general, the coefficients of the spectral series in terms of P 2 must 
vanish in order to fulfill (66a) - (66c). Since the diagonal coefficients b\ are equal to 
zero, the functions if^Tw) can be chosen arbitrarily. For the other functions, we find 



T A - - 

J 3.2 — 



«3 
b 2 . 



d J A 

dr. 



'4.2 



'5.2 



3 

Q3 «4 

Q3 «4 «5 

b 2 b 2 b 2 

3 4 5 



2.2 ; 
2 



d 1 JA 



M.3 



a4 



& 3 , 



~ 



4 



Li I w 

rl 2 ~ 

a JA J A 

d^.' ■■■■ ' hA 



a 5 



J_JA 

dr v 



(67a) 

(676) 
J 4 A 4 .(67c) 



Defining J a (t w ) ■= r Yn(T w )5 A y + S n (T w )5 A (f> with n > 2, where 7„(r t0 ) and 5„(tj„) 

are arbitrary functions of t w , we deduce from the recursive behavior of I A k(r w ) m 
(67a) - (67c) the general form of the time dependent coefficients of I a +1 (t w , y) as 

jA 



Jn+l( T w) 



I n+l.l( T w) = < 

where 

\n+k n 



n+1 



jn—l 



7 i T w ) 



for l = (n + l) 
for 2 < I < n , 



c\ ■= (-iy 



a TT fc b k 


ris=fe+i 


'2(s+2)(s-l)~ 
s+1 




3 fl s Ilt=3 


rit=fe+i 


t(t + 1) - k{k + 1) 



(68) 



(69) 



The last term in the right hand side of (69) can be expressed by factorials using the 
computer algebra program Maple and some properties)! of the Gamma function, 

2"- fc - 1 fc(n + 2)(2fc + 2)!(n-l)! 
' fc (fc + 2)!(n + fc + l)!(n-fc)! ' { ' 

We find the solution for U A with n > 1 by inserting I A into — 0, which 
leads to 



E (* + 2)(Z- 1)1 P,W(») , 



(71) 



where there is no summation performed over the index A on the right hand side, and 

ryA ._ „ X A , -l X A 0„++;„„ <E>(«) 
V 



sd A , + s" 1 ^. Setting TZ^' D) in (A.16) for n > 1 equal to zero and inserting 



I^+i an( i U%, gives rise to the solution (for n > 1) 

n+l 

*n+i(r w , y) = 



;=2 



J(Z + l)(i + 2)(i-l) 



n(n + 3) 



(72) 



In summary, the axially symmetric vacuum solution for with arbitrary n > 1 
is given by (3 n +i in (57), the coefficients I n +i.i in (68) for the power series (65), U A 
in (71), and $„+i in (72), respectively. For this solution each set of coefficients % is 
determined by two functions 7„+i(r, i ,) and S n +i(T w ), which can be prescribed freely 
along the geodesic c(r). 



t) Denoting with T(n) the Gamma function, we use r(n+l) = n! and r(n + 1/2) = 2- 2n Tr 1 / 2 (2n)\/n\. 
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4. Summary and discussion 

We studied the boundary conditions of the Bondi-Sachs metric functions & e 
{7, 5, (3, U A 1 $} at the vertices of null cones assuming that these vertices are traced 
by a time-like geodesic c(r) and that the observer emitting the null rays from this 
geodesic is an inertial observer. General requirements of these boundary conditions 
were found in three major calculational steps after assuming c(r) was contained in 
a convex normal neighborhood: In the first step, we constructed a metric in Fermi 
normal coordinates along c(t). In the second step, we defined affine null coordinates 
at c(r), where the radial coordinate is an affine parameter, and transformed the metric 
from Fermi normal coordinates to an affine null metric, which agrees with that of Ellis 
et al. (1985), Poisson (2004, 2005) and Poisson et al. (2006a, 2006b, 2010, 2011). In 
the third step, we calculated a Bondi-Sachs metric along c(r), by changing in the 
affine null metric the radial coordinate over to an areal distance coordinate r that 
defines the radial coordinate in the Bondi-Sachs metric. The boundary conditions 
of & at the vertex were found as Taylor series in terms of r with respect to a flat 
metric along c(r) where the expansion coefficients of the series are called correction 
coefficients 

Regularity at the vertex implied the following five general requirements on the 
power series of 

(i) The power series of & must start at r = with a certain positive power of r, in 
general 7, S, [3, $ are of 0{r 2 ) and U A of 0(r).ff . 

(ii) The coefficients ^„ have a rigid angular structure that is given by polynomials of 
harmonic and regular functions on 2-spheres centered on the vertices. 

(Hi) The polynomial coefficients of the harmonic base functions of ^ n carry strict 
numerical factors depending on the physical problem under consideration. 

(iv) The correction coefficient c € n , (n > 1) depends on the time derivatives of order 
(n — k) of the correction coefficient Ck, (1 < k < n) and these time derivatives 
must be finite. 

(v) Non-linear coupling occurs at lowest order near the vertex in the third-order 
correction coefficient of i.e. at 0(r 4 ) for 7, 5, (3, $, and at 0(r 3 ) for U A . 

Requirements (i) — (v) on the boundary conditions of & are kinematical regularity 
conditions, as they are derived from a regular metric along the geodesic tracing the 
vertices. The Einstein equations are not used in this establishing these conditions. 
The requirement that the ^ n have to obey Einsteins equations, this is a dynamical 
condition. 

We summarized these general requirements on the boundary conditions of J^" in 
a Vertex Lemma, which also answers question (I) from the introduction: 

Let the Bondi-Sachs metric functions, & ' , be represented by power series expansions 
with respect to the areal distance coordinate r at the vertex of a null cone on a time-like 
geodesic. If the coefficients of this radial expansion obey 

fflf the curve tracing the vertices is no geodesic but a general time-like curve, with a tangent 
vector u a , the lowest order terms depend differently on r, and one has to consider the acceleration, 
a M := u a \7 a u^, of the curve. Poisson et al. (2011) consider an affine null metric along such a general 
time-like curve. Transforming this metric to a Bondi-Sachs metric tells one how property (j) should 
be changed in that case. However, with our assumption of axial symmetry, it is natural to assume 
that the curve is a geodesic on the axis of symmetry. 
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(1) the general properties (i) — (v), 

(2) and the Einstein equations, 

then this power series expansion of & can be used to formulate boundary conditions 
for the metric functions & at the vertex of the null cone, and the vertex is a regular 
point in the null cone. 

We investigated the implication of the vertex lemma for axially axisymmetric vacuum 
space-times, which are the simplest space-times that allow us to demonstrate its 
relevant features. We also corrected and generalized the boundary conditions for 
such space times existing in the literature (Gomez et al. , 1994, Siebel, et al. , 2001). 
Since axially symmetric space-times contain a designated class of observers defined by 
time-like geodesies on the axis of symmetry, we introduce a Fermi normal coordinate 
system with respect to one of these axial geodesies c(r). 

We then proposed two approaches to answer question (II) of the introduction, 
i. e. how one can determine the boundary conditions of the Bondi-Sachs metric at the 
vertex. 

In the first one, we calculated the Riemann normal components in the Fermi 
normal coordinate system, and contracted these components with a null vector at 
c(t), its angular derivatives, and the tangent vector of c(r). These contracted 
Riemann normal component provided first-order correction coefficients of the Bondi- 
Sachs metric with respect to a fiat metric at the vertex. Although this approach works 
well for the first-order corrections, it is unsuited to find the higher order correction 
coefficients, because one would have to perform tedious full contractions of covariant 
derivatives of the Riemann normal components, and more importantly, because the 
approach itself cannot easily be extended to higher order. 

In the second approach we proposed the usage of a power series expansion of the 
Bondi-Sachs variables where the functions are expanded with respect to positive 
powers r, the coefficients depending on the time at the vertex and the direction angle 
of the emanating null ray. The lowest powers of r in these series were chosen to agree 
with those given in property (i) of the vertex lemma. We then showed that the other 
properties of the vertex lemma (ii — v) follow when integrating the Einstein equations, 
considering only the ^l — , ^2—, and ^3— correction coefficients. 

Considering property (ii) of the vertex lemma, we solved the partial differential 
equations for the and that they can in principle depend on the associated Legendre 
polynomials of first and second kind, P™(y) and Q l m (y), respectively. Since the Q™(y) 
are singular at the poles, y — ±1, we rejected them as possible solution. This is in 
agreement with the 'near- roundness' condition of Choquet-Bruhat et al. (2010). The 
polynomials Pf l (y) are the harmonic and regular functions mentioned in property 
(ii) of the vertex lemma. We further note that the polynomials, Qf l (y), were not 
encountered in the first approach where the Bondi-Sachs metric was determined 
directly from the Fermi metric, because the regularity at the poles is already assured by 
choosing the coordinate frame to be regular, i.e. choosing a Fermi normal coordinate 
system along c(r). 

We also confirmed properties (Hi) and (iv) of the vertex lemma, by exploiting the 
fact that the spectral expansion of n > 1, contains time derivatives of the indicated 
hierarchy (see table 2) and carries numerical factors following from the solution of the 
Einstein equations. 

Finally, we proved property (v) of the vertex lemma by showing that the solution 
for the ^3 coefficient depends linearly on the square of ^\ . 
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Table 2. The hierarchical dependence of the time derivatives of the free functions 
arising in the solution of the correction coefficients C2, and V, where the time 



derivatives with respect to t w 


are denoted by ovcrdots. 


correction free 


time 


coefficient functions 


derivatives 


72 (t w ), 52(t w ) 





c g 2 73 (t w ), 8 3 (t w ) 


72 ( r ™)> S 2 (t w ) 


% 74 (t w ), S 4 (t w ) 


73( T »): 72( r »). 63(t w ), 



For an axially symmetric vacuum space-time the solutions for each of the ^y — , 
^2 — , and coefficients depend on two time-dependent functions along c(t) which 
can be prescribed freely, i.e. both functions are not constrained by the vacuum 
Einstein equations. If these two free functions are (3— n) times differentiable with finite 
derivatives along c(r), the solutions of ^i, ^2 and ^3 provide the explicit boundary 
conditions for the Bondi-Sachs metric with a regular vertex along c(r) and for an 
axially symmetric vacuum. These six free functions correspond to the I = 2,3,4- 
multipoles of the electric and magnetic part of the Weyl tensor evaluated at the vertex 
in agreement with the findings of Thornc(I980), Hartle & Thorne (1985) and Zhang 
(1986), who found the metric expansions of a vacuum space-time along a timeline 
geodesic in deDonder coordinates. The boundary conditions found in the second 
approach approximate a regular vertex up to third-order corrections in Fermi normal 
coordinates and incorporate the non-linear coupling due to curvature at lowest order. 
Thus, they provide the answer to question (III) raised in the introduction. 

As an example for such non-linear boundary conditions let us consider the 
following six functions 

72{t w ) = K , j 3 {t w ) = , j 4 {t w ) = ^rK 2 , S 2 (t w ) = 5 3 (t w ) = 5 4 (t w ) = , (73) 

35 

where K > 0. Inserting these functions into the radial expansion coefficients in the 
solution for the n e {1, 2, 3} corrections in section 3.2.1, and comparing the 
thus obtained expressions with the radial expansion of the static, axially-symmetric 
vacuum solution, so-called SIMPLE (see Bicak et al. 1983, Gomez, et al. 1994 or 
Appendix B), in (B.7) - (B.12) shows that they agree. Hence, putting the functions 
(73) into our general solution for 'rfi, ^ 2 and ^3 approximates SIMPLE in a sufficiently 
small neighborhood of the vertex. Moreover, using the Bondi-Sachs metric at the 
vertex determined by our first approach, we can identify K with the Riemann normal 
component 6 x 7^1313, which gives rise to the (I = 2)— multipole of the electric part 
of the Weyl tensor. 

By linearizing the vacuum Einstein equations with respect to the coefficients 
^71, 1 < n < N, of we found how the hierarchical pattern of the time derivatives 
can be determined for arbitrary N, which provided us with the complete boundary 
conditions at a regular vertex for the Bondi-Sachs metric in linearized gravity and for 
axially symmetric vacuum space-times. The solution depends on 2A^ free functions, 7„ 
and 5 n , along the geodesic c(r). Regularity of the vertices along c(r) is guaranteed, 
if the two free functions in the coefficient ^„ are at least (N — n) times differentiable 
with finite derivatives along c(r). These 2N functions also determine the complete 
axially, symmetric initial data on the null cone for the linearized Bondi-Sachs metric 
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in vacuo. We claim that these 2N functions also determine the complete non-linear 
initial data on the cone to arbitrary order of N. Up to N = 3, this claim is backed up 
by our calculations. For N > 3, we learn from the coupling of the Riemann normal 
components in the Taylor series of normal coordinates at a space-time point (Schouten, 
1954, Thomas, 1991) that the non-linear coupling in a higher order coefficient occurs 
only between Riemann normal components and covariant derivatives determining 
lower order coefficients. Hence, our claim yields the answer to question (IV) raised in 
the introduction. We formulate it as a Conjecture: 

Suppose (^# 4 , g) is an axially symmetric vacuum space-time, c(t) a time-like 
geodesic on the axis of symmetry , JV W is a future (w = 1) or past (w = —1) 
null cone with its vertex on c(r), g(x a ) is a Bondi-Sachs metric on jY w 
with respect to coordinates x a — (t w , r, y, <f>), and its six metric functions 
J? are represented by the finite power series (28a)-(28e) in terms of the 
areal distance r with coefficients 1 < n < N evaluated at 6 '. Then the 
initial data on ,jV w for J? are given by 2N functions along c(r), which can 
be chosen freely. Regularity at requires that ^„ must be decomposed by a 
spectral series of associated Legendre polynomial P™+i{y) w ith respect to the 
angular variable y and that the two free functions in the coefficient ^„ are at 
least (N — n) times differ entiable and the derivatives of these functions are 
finite along c(r). 

The main consequence of the conjecture concerns the arbitrariness of the vacuum 
initial data that can be imposed on a null cone in Bondi-Sachs coordinates. If the null 
cone is assumed to have a vertex and if this vertex is assumed to be a regular point 
in the null cone, then one cannot impose any data. Instead the data have to obey a 
regular angular structure as it is provided by the solution of the Einstein equations 
for the coefficients of a Taylor series expansion of the metric at the vertex. The only 
freedom one has left in the choice of the data are free functions on the geodesic tracing 
the vertex. These free functions determine the highest I— multipole of the Legendre 
basis P ; m in each power of r in the Taylor series of the initial data at the vertex. They 
correspond to the /— multipolcs of the axisymmetric parts of the electric and magnetic 
parts of the Weyl tensor evaluated along the time-like geodesic. 

Our study raises further questions which should be investigated such as: How 
many free functions determine the initial data on a null cone with a regular vertex in 
space-times with no symmetry? What are the properties these free functions have to 
obey so that the initial data on the null cone describe an asymptotically flat space- 
time? What is the structure of the Bondi-Sachs metric at the vertex, if matter is 
assumed at the vertex and in its neighborhood? We plan to address the last question 
in an future publication. 
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Appendix A. Tensor components for the Einstein equations 

In the following, we give the components of lZ a fj that are used in Sect. 3 to calculate 
the solution of the vacuum Einstein equations up to third-order corrections and for 
the solution of the n th — order correction coefficients of linearized vacuum Einstein 
equations. We also introduce the following notations 

(A.1) 

s 2 e- 2 ~ t K yy + s^c^K^A cosh" 1 (26) , 



Tl(2D) ■= g AB TlAB , 



n 



(7) 



K 



(S) ■ = 



S G "^^yy & 6 "T^-(j^ 



sinh(2<5) - 2 cosh(25)^ . 



(A.2) 
(A.3) 



Appendix A.l. Ricci Tensor Components for the Third- Order Correction Coefficients 
The Ricci tensor components for the supplementary conditions are 



s 2 ($ 2 + 2/3 2 ),2, +6($ 2 + 2/3 2 ) + <^ s 2 ($3 + 2/3 3 ) 



12($ 3 + 2/3 3 ) 



Al32,T m +2§2,T W -Ul y \wr+\ s 2 ($ 4 + 2/3 4 ), y + 20($ 4 + 2/3 4 



2$ 3 ,r„ - 6/V„ - Ul yTa + (2$ 2 + 4/3 2 - 2 72 ) S 2 (<5 2 + 2h\ v 



+ stilly - 2($ 2 ,„ + 2^2,^)72,, + 2($ 2 ,,) 2 + 6(/3 2 ,,) 2 + 6$ 2 , y /3 2 ,, 



-\(ut) 2 



s 2 - 



UfU% y + -(t/^) 2 ] S 4 - ($ 2 , v + Ut yy + 8fc, v )U? 



s - 2 uf 



r- x s 2 ll 



+ (96^ 2 - 2Ul y )<5>2 + 32 [($ 2 ) 2 + 2(/? 2 ) 2 
- (Ul y ) 2 2s-VK) 2 jr 2 + 0(r 3 ) 

2(02 + <S>2),yS 2 + 3Uf] W + { [3(03 + *3) - fo,T v ] + 6C/| 



(A.4) 



•s 2 + (s 2 73,T„) l!y 



+ (S 2 7 2,rj,, - \ut, Tw }r + I [4($ 4 + &) + for* 

+ {«*h, v ),vV\ + 4 [( 2 *2 + 3/3 2 - 72 )($2 + &),,, + 4(3C/f - y^ y )<y 2 ] s 2 

+ [j 2 [^ w - 2C/ft/^] s 4 + (l0$ 2 + 2f3 2 + 14 72 - y U Iv) U " + 10C/ 3 

~ Ulr w ~ 4ys' 2 (C/|') 2 }w;r 2 + 0(r 3 ) (A.5) 



(s 4 ^ + 2Uf] w + { [^ 4 C/t),, + (s 2 S 2 ,rJ 



r + 



~(s 4 Ul y ), v + (s 2 5 3 , T J, y - (s 4 7 2), y U? 2(s i uiy), vl 2 



+ 9U+ - Ul Tw - 2yUfl3 2 , y - 2 S 2 (7^ 72 ,, + fa - \u\ tV - 10 72 + 10$ 2 ) uf 
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+ 2S 2 , y Ul y }+0(r 3 ) 

The Ricci tensor components for the hyper-surface equations are 

16/3 4 - 8( 72 2 + <5 2 2 )lr 2 +0(r 3 ) , 
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2yUf + 75 2 ] s - 2 U? 



(A.6) 



n 



r^ x s 2 Tl r 



= 8/3 2 + I2(3 3 wr + 



2Uf + 2{ S 2 l2 ). y + 5C/| - s 2 (3 3 , y + 3( S 2 73 ),y 



(A.7) 



rw 



9U$ - 2s 2 /3 4 , y - 4 S 2 7272 ^ + 4(^72)^ + (6U+ - 4S 2ty )s 2 5 2 
+ 6( l2 -(3 2 )Uf]r 2 +0(r 3 ) , (A.8) 

r _1 s -2 ^ = 2Uf + 2s- 2 {s 2 5 2 ) ty + huf + 3s~ 2 (s 2 5 2 ) J rw' 1 

+ \9U$ + 4 S - 2 (s 2 S 4 ), y + 6s- 2 UfS 2 - 6( 72 + fo)U}\ r 2 + 0(r 3 ) , (A.9) 
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(215) 



= - 12$ 2 + 5Ul y - 2(/3 2 , yS 2 ) y - 8/3 2 + 2 S - 2 ( S 4 72,,) >y - 4 7 2 
+ [ - 16$ 3 + 6Ul y - 2(p 3iV s 2 ) y - 12/3 3 + 2 S - 2 ( S 4 7 3,y)^ - 4 7 3 

+ | - 20$ 4 + 7Ul y - 2{^ y s 2 ) y - 16ft + 2 S - 2 ( S 4 74 ^) _„ - 4 7 4 
- 4 S 2 ( 7 2 + 51) 20$ 2 -\(^) 2 - l(sUf) 2 4 S - 2 ( S %, y )6 2 



2S 2 (f3 2 ,y - 2 l2 ,y)(3 2 ,y + 4 (S 2 (3 2 ,y),y - S " 2 {.^2 ,y ) , y 



72 



+ 4(72 + Si) - 2(8$ 2 + 5C/^)/3 2 \r 2 + 0{r 3 ) 



(A.10) 



and for the evolution equations 

r- 2 K {l) = - 12 72 - 2s 2 ^ yy + 2,s 2 (s- 2 U\ ),„ + [l2 72 , rra - 24 73 - 2s 2 /3 3 , w 

+ 4 S 2 (s- 2 C/|), y ] wr + |l673, r „ - 40 74 + 5s 2 ( S - 2 [/|), 2/ - 2s 2 ^ yy 

1 fUf\ 2 

- ^2I72Ti*li 1 72,j- 

1 



- e* 2 /^- 2 !/?),, + 4([^ iW - 10$2 - 4/? 2 ) 72 + 14C/f72,, - \(^-) 



4 72 /3 2 , w - 2^ iV + U6 2 U« y + -(Uf) 2 \ S 2 ^r 2 +0(r 3 ) , (A.ll) 
r- 2 K (s) = 12S 2 - 3s 2 Ul y + ( - 126 2 , Tw + 24S 3 - 4s 2 uQwr 

+ | - 16<y 3 ,r„ + 40<5 4 - 5s 2 ut y + [(6& + 1472)^ - 4/3 2 , TO 5 2 ] s 2 

+ (I/? - 14J 2 , W )I/* + 4(4/3 2 - + 10$ 2 )l<5 2 }r 2 + 0(r 3 ) . (A.12) 
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The Ricci quantities for the main equations that are linearized with respect to are 
assumed to be of the form 

N N 



n 



n=0 n=0 

where 7£(.) G {^(7)) T^(2D)}- The relevant non-zero coefficients for the 

surface equations read 



ft 



(n) 

(2D) ~ 



{An + 8)/3„ +2 

Jn(n + 3)f7« + (1 - n)/?„+i. y + (n + 1) f s 2 7 „ +1 ) 
^ v ' ,y 

X -n(n + 3)17* + S - 2 (n + 1) (s 2 5 n+ i) 

-4(n + 3)$„ +2 - 4(n + 2)/}„+2 + (n + 5)l^ +1 ^ 

- 2[s 2 l3 n+ 2.y ) + 2s~ 2 (s 4 7„ +2 . y J -47„ +2 
v 7 ,1/ v y ,y 



n > , 
n > 1 , 

n > 1 , 



hyper- 

(A.13) 
(A.14) 

(A.15) 



and the evolution equations for 7 and (5 are given by 



n 



K 



(2) 
(7) 

(2) 
(-5) 



-12 7 2 - 2s 2 /3 2 , : 
125 2 - 3t?Ut iy 



3s 



! («- a t^) 



n > , (A.16) 

(A.17) 
(A.18) 



and by 



nf+ 1] = 2(n + 1) [27 n , TM - (n + 2) 7 „+i] + s 2 (n + 2)(V 2 t/^) ^ - 2s 2 /3„ +1 

ft[™ +1) - 2(n + 1) [(n + 2)5„+i - 2<5„, r J - S 2 (n + 2)17* „ . 
for n > 1, respectively. 



;/;/ 



(A.19) 
(A.20) 



Appendix B. The vacuum solution SIMPLE 

The static vacuum solution SIMPLE is the only known explicit non-linear and axially 
symmetric solution of the Einstein equations in Bondi coordinates (Bicak et al. 1983). 
SIMPLE is boost and rotation symmetric (Bicak et al. 1984), and depends on one 
free parameter a. The solution given in Gomez et al. (1994) reads in our coordinates 
x a = (t w , r, y, (f>) and nomenclature of variables 



S(o, r, y) = + a 2 r 2 s 2 (y) 
7(r,y) = ln(l + £(a,r,y))-ln(2) 



P(r, y) 
U y (r, y) 

$(r, y) 
S(r, y) 



In 

a 2 rys 2 (y) 
S(o, r, y) 
1 



1 + E(a, r, y) - - In £(a, r, y) - In 2 , 



In 



2aW(y)-aV + l 



In 



1 + S(o, r, y) 



+ ln2 



uHr, y) = . 



(B.l) 
(B.2) 

(B.3) 

(B.4) 

(B.5) 
(B.6) 
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Defining a = (12K) 1 / 2 and expanding the Bondi functions of SIMPLE near r = 
gives 



7(r,y) =r 2 KPi{y) + 
5(r,y) =0 , 
P(r, y) 



2 lK 2 Pi{y) + ^K 2 Pl{y) 



9 -K 2 P 2 (y)-^K 2 Pt(y) 



UV(r, y) = -ArKs{y)Pl{y) + 
$(r, y) = -6r 2 ^P 2 °(y) - 



96 



r 4 + 0(r 6 ) , 
144 



^ + ^P 2 0( y )-^pO( y) 



r 4 + G(r 6 ) , (B.7) 
(B.8) 
(B.9) 

S ( y y + 0(r 5 ), (B.10) 
(B.ll) 
r 4 + 0(r 6 ) , 

(B.12) 



where we expressed powers of y in terms of associated Legendre polynomials. 
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